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Recent experiments by Hayashi et al. [Phys. Rev. Lett. 91, 226804 (2003)] demonstrate coherent 
oscillations of a charge quantum bit (qubit) in laterally defined quantum dots. We study the intrinsic 
electron-phonon decoherence and gate performance for the next step: a system of two coupled charge 
qubits. The effective decoherence model contains properties of local as well as collective decoherence. 
Decoherence channels can be classified by their multipole moments, which leads to different low- 
energy spectra, ft is shown that due to the super-Ohmic spectrum, the gate quality is limited by 
the single-qubit Hadamard gates. It can be significantly improved, by using double-dots with weak 
tunnel coupling. 
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I. INTRODUCTION 

In recent years, the experimental progress in analyz- 
ing transport properties in double quantum dotai has 
lead to the fabrication of double dot structures with only 
one electron in the whole systeroSi^. This well-defined 
situation permits, although it is strictly speaking not 
necessaryii, to use quantum dot systems as quantum bits 
(qubits). In order to define qubits in lateral quantum 
dot (QD) structures, the two degrees of freedom, spin 
and charge, are naturally used. For spin qubits^, the in- 
formation is encoded in the spin of a single electron in 
one quantum dot, whereas for the charge qubili^iii^ the 
position of a single electron in a double dot system de- 
fines the logical states. Similar ideas can also be applied 
to charge states in Silicon donors^. Both realizations are 
interconnected: interaction and read-out^ of spin qubits 
are envisioned^ to be all-electrical and to make use of the 
charge degree of freedom. 

Although the promises of spin coherence in theorjiiS 
and in bulk measurementsii are tremendous in the long 
run, it was the good accessibility of the charge degrees of 
freedom which lead to a recent break-through^, namely 
the demonstration of coherent oscillations in a quantum 
dot charge qubit. In this experiment, three relevant de- 
coherence mechanisms for these charge qubits have been 
pointed out: a cotunneling contribution, the electron- 
phonon coupling, and l//-noise or charge noise in the 
heterostructure defining the dots. 

Recent theoretical results-^ predict that the cotunnel- 
ing contribution can be very small, provided that the cou- 
pling between the dots and the connected leads is small. 
Thus, cotunneling is not a fundamental limitation. This, 
however, means that initialization and measurement pro- 
tocols different from those of Ref. Q are favorable^. 

Other theoretical worksi^ti^ii^iiSiii already describe 
the electron-phonon interaction for a single charge qubit 
in a GaAs/AlGaAs heterostructure. Moreover, also elec- 
tronic Nyquist noise in the gate voltages affects the qubit 



systen>W-. Note that the physics of the electron-phonon 
coupling is different and less limiting in the unpolar ma- 
terial Si^^, where the piezo-electric interaction is absent. 



II. MODEL 

In this article, we analyze the decoherence due to the 
electron-phonon coupling in GaAs, which is generally as- 
sumed to be the dominant decoherence mechanism in a 
coupled quantum-dot setting. The recent experimental 
analysis shows that the temperature dependence of the 
dephasing rate in the experimen1?i can be modeled with 
the Spin-Boson model and hence is compatible with this 



assumptio] 
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We develop a model along the lines of 



Brandes et almii^ to describe the piezo-electric interac- 
tion between electrons and phonons in lateral quantum 
dots. Thereby, we assume the distance between the two 
dots to be sufficiently large and the tunnel coupling A to 
be relatively small, which is a prerequisite for the validity 
of the model. The Hamiltonian for a system of two dou- 
ble dots with a tunnel-coupling within the double dots 
and electrostatic coupling between them, see Fig. ^ can 
be expressed as2S 



^total — ^^sys + -ffbath + ^^int , 



(1) 



where 



^sys = ~ X! ^ {£i^z,i + Ai(Ta;^i) - fc(T^,i (g) a^^2 (2) 



bath 



(3) 



refer to the qubits and the heat bath, respectively, q is 
the phonon wave number. The system-bath interaction 
Hamiltonian Hi^t depends on details of the setup such as 
the crystalline structure of the host semiconductor and 
the dot wave functions. We will distinguish between the 



two extreme cases of long correlation length phonons re- 
sulting in coupling of both qubits to a single phonon bath, 
or two distinct phonon baths for short phonon correla- 
tion length. The former case is more likelyS^ and ap- 
plies to crystals which can be regarded as perfect and 
linear over the size of the sample, whereas the latter 
case describes systems that are vstrained or disordered 
and double quantum dots in large geometrical separa- 
tion. The correlation length has to be distinguished from 
the wave length: The former indicates, over which dis- 
tances the phase of the phonon wave is maintained, i.e., 
over which distance the description as a genuine stand- 
ing wave applies at all, whereas the latter indicates the 
internal length scale of the wave. 



A. One common phonon bath 

In the case of a single phononic bath with a very long 
correlation length coupling to both charge qubits, -ffint 
can be written as 
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(4) 



The coefficients aq^i and [3q^i describe the coupling of a 
localized electron (one in each of the two double dot sys- 
tems) to the phonon modes. They are given by 



aq,i = \q{l,i\e'^^\l,i), 
Pq,, = A,(r,z|e^«->,z), 



(5) 
(6) 



where the \l,i) and \r,i) denote the wavefunctions of 
the electrons in the left or right dot of qubit i.We as- 
sume these wavefunctions to be two-dimensional Gaus- 
sians centered at the center of the dot, as sketched in 
Figure n These states approximate the ground state in 
the case of a parabolic potential and small overlap be- 
tween the wavefunctions in adjacent dots. The coefficient 
Xq is derived from the crystal propertied. 

Henceforth, we investigate the case of two identical 
qubits. Due to the fact that the relevant distances are 
arranged along the a;-direction, we obtain the coupling 
coefficients 



aq.i 

aq,2 



^ y^^^iq(l/2+d)^-q^a^/i ^ 



(7) 
(8) 
(9) 

(10) 



Here, q is the absolute value of the wavevector q. The 
second exponential function in each line is the overlap 
between the two Gaussian wavefunctions. 

This two-qubit bath coupling Hamiltonian is quite 
remarkable, as it does not fall into the two standard 
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FIG. 1: (Colour online) Sketch of the two coupled identical 
charge qubits realized in a lateral quantum dot structure, d = 
100 nm is the distance of the dot centers in one qubit, I — 
200 nm is the distance between the right dot center of qubit 1 
and the left dot center of qubit 2. The width of the Gaussian 
wavefunction of an electron in each dot is cr = 5 nm. The 
values chosen for the distances d and / are slightly smaller 
than in experimental realizations^*- in order to provide a lower 
bound for the decoherence times. In principle, there could be 
tunneling processes between both qubits, i.e., the QDs two 
and three in the chain, but we assume that the coupling is 
pinched off by applying appropriate gate voltages. The gray 
box between the qubits indicates that there is no tunneling 
between the qubits. 



categorie s usually treated in literature (see, 
'24"'25'26] and references therein): On the 



e.g., Refs. 
one hand, 

there is clearly only one bath and each qubit couples 
to the bath modes with matrix elements of the same 
modulus, so the noise between the qubits is fully corre- 
lated. On the other hand, the Hamiltonian does not obey 
the familiar factorizing collective noise form -ffsB.coU = 
^system ® ^bath- Such a form would lead to a high de- 
gree of symmetry and thus protection from the noise 
couplingSiS^, however, the Hamiltonian Hint, eq. 01 can- 
not be factorized in such a bilinear form. It is hence 
intriguing to explore where in between these cases the 
physics ends up to be. This is in particular important 
for finally finding strategies to protect the qubits against 
decoherence, and for estimating the scaling of decoher- 
ence in macroscopic quantum computers. 

In order to obtain the dynamics of the reduced density 
matrix p for the coupled qubits, i.e., for the degrees of 
freedom that remain after the environment is traced out, 
we apply Bloch-Redfield theorji2L22i2S. It starts out from 
the Liouville-von Neumann equation ihp = [H, ptot] for 
the total density operator. A perturbational treatment 
of the system-bath coupling Hamiltonian /fjnt results in 
the master equation 



P 



■[Hsys,p]- 



1 



dr trs[i7i„t, [Hinti-r), ps]], 

where ps = exp(— /3iJB)/Z denotes the equilibrium den- 
sity matrix of the bath. Evaluating the trace over all 
bath variables, tr^, and decomposing the reduced density 
operator into the eigenbasis of the unperturbed system 



3 



Hamiltonian, we obtain2Si2£ 

where ujnm = En — Em ■ The first term on the right hand 
side describes the unitary evolution and the Redfield re- 
laxation tensor Rnmke incorporates the decoherence ef- 
fects. It is given by 

nnmki — Oim 2_^^ nrrk~^ "nk 2_^^ Irrm ^ imnk ^ imnk^ 
r r 

(13) 

where the rates T^^^ are determined by Golden Rule 
expressions^Mi, see Eqs. (jSOJ and |2lJ, below. The 
Redfield tensor and the time evolution of the reduced 
density matrix are evaluated numerically to determine 
the decoherence properties of the system due to a weak 
electron-phonon coupling. Note that in addition, Ohmic 
electronic noise can be taken into account by employ- 
ing the spectral function^^ ^s(w) = ^ohmic('-^) + J{'^), 
where J{lij) contains only the phonon contribution. It 
is also possible to take l//-noise in the quantum dot 
system into account in the same way. The l//-noise 
essentially determines the magnitude of the dephasing 
part of the decoherence. Thus, it is in turn possible to 
impose for the zero frequency component J(0) the ex- 
perimental value of the dephasing rates or a value from 
a microscopic model22.. However, in many cases it turns 
out to be non-Markovian and/or non-Gaussian, leading 
to non-exponential decay, which can neither be described 
by Bloch-Redfield theory nor parameterized by a single 
rate. 



In order to compute the rates, the electron-phonon in- 
teraction Hamiltonian has first to be taken from the lo- 
calized representation to the computational basis, which 
is straightforward. To compute Bloch-Redfield rates, it 
is necessary to rotate into the eigenbasis of the system. 
After this basis change, the spectral densities Jemnk{^) 
are calculated along the lines of Ref. j22] as 

Jimnkicj) - {{B-^CB)^^ {B-^CB)J„ (14) 

where B is the matrix for the basis transformation from 
the computational basis {|00), |01), |10), |11)} to the 
eigenbasis of the system and {■)q denotes an averaging 
over all phonon modes q with frequency uj. The matrix 
C is diagonal in the computational basis, C = diag(aq.i — 

Pq.l + ag,2 - /?g,2, - Pq.i + aq,2 - /3q,2, ^g,! " (3q,l + 

Olq,2 - /3q,2, Ckga - (3q^l + Q;g,2 - f3q,2)- 

The explicit derivation shows that it is most convenient 
to split the total spectral function Jimnk{^) [see Eq. (|14|l ] 
into odd and even components 

Jo{oj) , (15) 

where the prefactors eemnk and oimnk of the even/odd 
part of the spectral function are matrix elements coming 
from the basis change from the computational basis to 
the eigenbasis of the system and 

Je/o{(^) = J^\aq,l-l3q,l±aq,2Tf3q,2fS{uJ-UJq). (16) 

q 

They evaluate to 



Trhujg 



2 - 2 — sm — =F — sm — ± 2 — — sm =F — - — sm 

UJ \L^d J UJ \UJl J UJ \UJl+dJ UJ 



yUJl+2d, 



(17) 



where g = 0.05 is the dimensionless electron-phonon cou- 
pling strength for the commonly used material GaAsSLS^ 
and cs the speed of sound. The different frequen- 
cies represent the distances in the system: uJd = Cs/d, 
uji = Cs/l, ujd+i = Cs/{d + l) and uj2d+i = Cs/{2d + l), 
and ujc = Csjcr. This structure can be understood as fol- 
lows: The electron-phonon interaction averages out if the 
phonons are rapidly oscillating within a dot, i.e. if the 
wavelength is much shorter than the dot size — this pro- 
vides the high-frequency cutoff at uJc- On the other hand, 
long-wavelength phonons do not contribute to decoher- 
ence between dots i and j, if the wavelength is much 
longer than their separation because then, the energy 
shift induced by the phonon displacement will only lead 
to a global phase. Furthermore, we can approximate the 



leading order at low frequencies as 

Uuj) = ^-^^^ + 0{uj% (18) 
3cf 

JoH = j^^4 -io^ + 0{uj^). (19) 

This different power-laws uj^ to uj^ can be understood 
physically as illustrated in Figure |21 "Even" terms are 
the natural extension of the one-qubit electron-phonon 
coupling, adding up coherently between the two dots. 
In the "odd" channel, the energy offset induced in one 
qubit is, for long wavelengths, cancelled by the offset in- 
duced in the other qubit. Thus, shorter wavelenghts are 
required for finding a remaining net effect. An alterna- 
tive point of view is the following: The distribution of 
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even contribution 




FIG. 2: Illustration of the even (top) and odd (bottom) con- 
tributions to the total rates. Filled circles indicate occupied 
dots. For long-wavelength modes, the energy shifts induced 
by underlying phonons in the two dots add up coherently in 
the even case but cancel in the odd case. Note, that moving 
charges from the black to the white dots changes the dipole 
moment in the even but not in the odd case. 
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For (jjij 0, these rates vanish due to the super-Ohmic 
form of the bath spectral function. From this, we find the 
time evolution of the coupled qubit system and finally 
also the gate quality factors. 



B. Two distinct phonon baths 

When each qubit is coupled to its own phononic bath, 
the part of the Hamiltonian that describes the interaction 
with the environment Hi^t is given by 

Hint = + + (oiqi - M^^A 

91 

X(4l 12 + 

92 

x(4+c_,j0ii. (22) 

This scenario can be realized in different ways: One can 
split the crystal into two pieces by an etched trench. Al- 
ternatively, if there is lattice disorder and/or strong non- 
linear effects, the phonons between the dots may become 
uncorrelated. 

The calculation of the coupling coefficients works in a 
similar way, but there are two different indices qi and q2 
to represent the phononic baths of each qubit 



FIG. 3: (Colour online) Spectral functions Je,o{i^) in the case 
of one common phonon bath for the fixed parameters Cs — 
5000 m/s, g = 0.05, d = 100 nm, I = 200 nm and a = 5 nm. 
Inset: zoom for small frequencies. 



the two charges can be parameterized by a dipole and 
a quadrupole moment. The "even" channel couples to 
the dipole moment of the charge configuration similar to 
the one-qubit case. The "odd" channel couples to the 
quadrupole moment alone (see Figure Thus, it re- 
quires shorter wavelengths and consequently is strongly 
supressed at low frequencies. This explains the different 
low-frequency behavior illustrated for realistic parame- 
ters in Figure |3| Thus, we can conclude that for small 
frequencies the odd processes are suppressed by symme- 
try — even beyond the single-dot supression and the sup- 
pression of asymmetric processes. 

With these expressions for the spectral densities, one 
can proceed as in Ref. 26] and determine the rates that 
constitute the Redfield tensor to read 

J imnk iy^nk ) 
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(23) 

(24) 
(25) 
(26) 



The expression for the spectral functions Jtmnkiy) 
turns out to be exactly the same as the one in the last 
section with the only difference that instead of a^.i, the 
coupling between electrons and phonons is now expressed 
as Qig. (with i = 1,2 for both qubits). Therefore, in or- 
der to obtain the spectral density Jenmki^^), one has to 
average over two distinct baths, i.e. 



J Imnki^^) - {{B ^CB)^^[B ^CB)^Jq^^q^. 



(27) 



Again, we find two different functions that we name in 
the same way as in the previous section, Je(ci^) and Jo{lu), 
which are given by 
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FIG. 4: (Colour online) Spectral functions Je,o(a;) in the case 
of two distinct phonon baths for the fixed parameters Cs = 
5000 m/s, g = 0.05, d = 100 nm, / = 200 nm and cr = 5 nm. 
Inset: magnification for small frequencies. 



The prefactors from the basis change also enter the ex- 
pressions for the rates in the same way as in the last 
section. The spectral functions Je^o(w) are plotted in 
Figure^ the inset depicts the proportionality to w'^ for 
small frequencies. 



III. GOLDEN RULE RATES 

We proceed as in Ref. and determine the Golden 
rule rates that govern the Redfield tensor. Thereby, we 
find both the time evolution of the coupled system and 
the gate quality factors. 

Let us first discuss the impact of this particular bath 
coupling on the dephasing and relaxation rates. The de- 
coherence rates, i.e., the relaxation and dephasing rates, 
are defined according to Tji = "X^n^m where A„ are 
the eigenvalues of the matrix composed of the elements 

Rji .71 ^ra .m ^ Tl^TTL 1, ... ,4, and Tip^^ — ^^Rn.7n,n.7n 



> \R 



and in 



for non-degenerate levels \ll!„ 
the absence of Liouvillian degeneracy, |wrim — i^kil > 
\Ra,b,c,d\ a, b, c, d, e {k, I, m, n], respectivelji^i. 

As a reference point, we study the rates in the uncou- 
pled case. In this case, and in the absence of degeneracies 
between the qubits, there is a clear selection rule that the 
environment only leads to single-qubit processes, i.e., de- 
coherence can be treated at completely separate footing. 
As a result, all rates are identical between the qubits. To 
make this obvious, we rewrite the original Hamiltonian in 



^int - 



iq{l+d)/2 



(29) 



which — besides a phase factor which is meaningless for 
single-qubit transitions — is identical to the standard 
electron-phonon Hamiltonian for double dots22i. 

Figure |S1 shows the temperature dependence of the en- 
ergy relaxation rate F/j and the two dephasing rates F^^g 
and F^24 compared to the single qubit relaxation and de- 
phasing rates. In this notation, F^^ . is the rate at which a 
superposition of energy eigenstates i and j is decays into 
a classical mixture. We considered the following three 
cases, characterized by values on the matrix element rel- 
ative to a characteristic system energy scale Eg', (a) 
large difference of the Ei and (i = 1,2) between both 
qubits and no coupling between the qubits (ei = Ai = 
{l/40)Es, 62 = A2 = -(21/40)Ss and coupling energy 
K ^ 0), (b) small asymmetry between the parameters 
for both qubits and no coupling (ei = Ai = — (l/2)i?s, 
£2 = A2 = ~{21/A0)Es and K ^ 0), and (c) with- 
out asymmetry between the qubits and a rather strong 
couphng between the qubits (ei = Ai = —{l/2)Es, 
£2 = A2 = -{l/2)Es and K = lOEs). One generally 
would expect a different value of the distance between 
the dot centers in the qubits d, when the tunneling cou- 
pling is varied. However, in our case of the dot wave- 
functions which overlap only in their Gaussian tails, this 
effect is very small (below 1 nm for a change in the tun- 
neling amplitude A of approximately ~ {l/2)Es) for the 
lengthscales that we are considering. Note, that in Ref. 
a substantial change of A over more than an order of 
magnitude was obtained experimentally by a rather mild 
adjustment of the gate voltage, so it is consistent that a 
small change of A can be achieved by a tiny adjustment. 
Therefore the value d = 100 nm is used for the electron- 
phonon coupling encoded in and Jo in all cases. 

For case (a), we find that all rates are for all temper- 
atures larger than the single qubit rates, as one would 
expect33. In more detail, for the single bath case, the 
ratio of the relaxation rates is approximately 1.9, the ra- 
tio of the single-qubit dephasing rate and the two-qubit 
dephasing rate Tc/j^i is around 0.9 and for the dephasing 



rate F, 



013 I 



the ratio is 1.0. The behaviour of the even 



and odd parts of the spectral function in the single bath 
case can be explained from the spectral function Fig. 13 
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FIG. 5: (Colour online) Temperature dependence of the relaxation and dephasing rates normalized by the single-qubit relaxation 
and dephasing rates. The two-qubit relaxation rate is given by the trace of the relaxation part of the Redfield tensor in secular 
approximation. The energy scales for the two-qubit transitions 1 <-> 3 and 2 <-> 4 are comparable to the single qubit energy scale, 
the characteristic qubit energies are Es = (1/8) GHz. The different cases are (a) ei = Ai = (l/40)-Es, £2 = A2 — — (21/40)i5s, 
and coupling energy K = 0, (b) ei = Ai = -{l/2)Es, £2 = A2 = -(21/40)^^, and K = 0), and (c) £1 = Ai = -{l/2)Es, 
£2 = A2 = — (l/2)i5s and K = lOEa. Note that cases (a) and (b) model uncoupled qubits, especially for case (a) the 
overall relaxation rate for the two-qubit system is approximately twice the single-qubit relaxation rate when calculated for the 
dominating larger energy scale of the two-qubit system (£2 = A2 — —(2 1/40) £3). 



for small w one finds that Jo < Je- For the case of large 
frequencies, however, the even part of the spectral func- 
tionincreases and even dominate beyond the threshold 
uj ^ ujd- Overall, it is found that in the case of a single- 
bath the decoherence effects are significantly suppressed 
compared to the two-bath scenario. For the two-bath 
case, the ratios are for the relaxation rates approximately 
3.9, for the dephasing rate T^^^ around 1.9 and for the 
dephasing rate F^^g it is 2.0. Note that for the two-bath 
case Je < Jo always and for the case where both tun- 
nel matrix elements in the Hamiltonian vanish, the rate 
vanishes, too. 

After decreasing the asymmetry between the two 
qubits as in case (b), the rates decreased but are still 
comparable with the single qubit rates, besides the last 
dephasing rate F^^^. This can be understood, if one con- 
siders the energy spectrum of the eigenvalues of the sys- 
tem Hamiltonian. In cases (a) and (b) there is signifi- 
cant difference between the qubits, so it is straightfor- 
ward to map the two-qubit rates onto the corresponding 
single qubit rates and they are largely determined by 
single-qubit physics. In case (c), we consider a fully sym- 
metric case in the qubit parameters, but with a finite 
and large coupling between the qubits. This coupling 
lifts the degeneracy but makes the rate a generic two- 
qubit rate which belongs to a relatively robust transition 
with small transition matrix elements for the single bath 



case. At high temperatures, these symmetry-related ef- 
fects wash out as discussed in Ref. [s^. However, the 
high-temperature rates do not coincide with the single- 
qubit rates, as the underlying energy scales are still dif- 
ferent and in generally larger for the two-qubit situation. 

Overall, the ratio of the two-qubit and single-qubit re- 
laxation rates decreases for increasing temperature due 
to the reduction of correlation effects in the double dot 
system, besides case c), where a symmetry based on the 
underlying Hamiltonian becomes important. 



IV. QUANTUM GATE PERFORMANCE 

For the characterization of the quantum gate perfor- 
mance of this two-qubit system, it is necesssary to intro- 
duce suitable quantifiers. Commonly, one employs the 
four gate quality factors introduced in Ref. [33 ; fidelity 
J-', purity V, quantum degree Q, and entanglement capa- 
bility C to chararcterize a gate operation within a hostile 
enviroment. 

The fidelity, i.e., the overlap between the ideal prop- 
agator and the simulated time evolution including the 
decoherence effects, is defined as 



(*in|L/tpo,t?7|^'i„), (30) 
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where the bar indicates an average over a set of 36 unen- 
tangled input states l^'in) = IV'j) with i, j = 1, . . . , 6. 
The 6 single-qubit states ji/'i) are chosen such that they 
are symmetrically distributed over the Bloch sphere, 

|V'l) = |0), 1^2) = |^3,...,6) - ^ (31) 

where (j) — 0, 7r/2, tt, 37r/2. Here, U is the ideal unitary 
time evolution for the given gate, and pout is the reduced 
density matrix resulting from the simulated time evolu- 
tion. A perfect gate reaches a fidelity of unity. The purity 
V measures the strength of the decoherence effects, 

V = triple,). (32) 

Again, the bar indicates the ensemble average. A pure 
state returns unity and for a mixed state the purity can 
drop to a minimum given by the inverse of the dimension 
of the system Hilbert space, i.e. 1/4 in our case. 

If the density operator p describes an almost pure 
state, i.e., if the purity is always close to the ideal value 
1, it is possible to estimate the purity loss during the gat e 
operation from its decay rate along the lines of Ref. [33 ■ 
Thereby, one first evaluates the decay of (c?/dt)trp^ for 
an arbitrary pure qubit state p = From the basis- 

free version of the master equation ({TT|l , follows straight- 
forwardly 

— trp^ = --2 / dr tr5+B[i7i„t, [iJint(-r), p ® Pb]]p- 
at n Jq 

(33) 

By tracing out the bath variables, we obtain an expres- 
sion that contains only qubit operators and bath corre- 
lation functions. This depends on the state 1-0) via the 
density operator. Performing the ensemble average over 
all pure states as described in the Appendix 1X1 we obtain 

^^¥{WTT)lo '^^tr([i?i„t,i?int(-r)]+)Bxq, (34) 

where iV = 4 denotes the dimension of the system Hilbert 
space of the two qubits. We have used the fact that 
triJint = 0. Although the discrete and set of states em- 
ployed in the numerical computation is obviously differ- 
ent from the set of all pure states, we find that both 
ensembles provide essentially the same results for the pu- 
rity. 

If the bath couples to a good quantum number, i.e., 
for [iJgys, ^^int] = 0, the system operator contained in 
the interaction picture operator Hi-at{—T) remains time- 
independent. Then, the r- integration in (|34|l is effectively 
the Fourier transformation of the symmetrically ordered 
bath correlation function in the limit of zero frequency. 
Thus, we obtain 

V ^ ^ — lim V Ji(tj)coth4^, (35) 



where 

J^i'^) = \'^\°'q-i~ f^q-A'^^^^ ~ "^q) (36) 

denotes the spectral density of the coupling between 
qubit i and the heat bath(s). 

In the present case of a super-Ohmic bath, the limit 
— > results for the coupling to a good quantum num- 
ber in 7^ — 0. This means that whenever the tunnel 
coupling in the Hamiltonian (jSJ is switched off, i.e. for 
Ai = A2 = 0, the purity decay rate vanishes. Thus, 
we can conclude that the significant purity loss for the 
CNOT operation studied below [cf. Eq. H41() ]. stems from 
the Hadamard operation. This is remarkably different 
from cases with other bath spectra: For an ohmic bath, 
for which Ji{uj) oc uj, expresion (|35|) converges in the 
limit a; — > to a finite value. By contrast, for a sub- 
ohmic bath, this limit does not exist and, consequently, 
the purity decay cannot be estimated by its decay rate. 
During the stage of the Hadamard operation, A2 = A 
while Ai — 0. Then, the interaction picture versions of 
the qubit-bath coupling operators read 

c^zA-t) = o-^j, (37) 
o'z,2{-t) = 0-2^2 cos(Ar/^) - (7j^^2 sin(AT/?i). (38) 

In the case where both qubits couple to individual en- 
vironments, the expression for the change of the purity 
can be evaluated for each qubit separately. For qubit 
2, we still have a coupling to a good quantum number, 
while for qubit 1, the appearence of cos(Ar/?i) results in 
a Fourier integral evaluated at the frequency A/h. Thus, 
we finally obtain 

For the super-Ohmic bath under consideration [see eqs. 
((TH|l and ((T^ ]. the first term in Eqn. vanishes. 

In the case of one common heat bath, the estimate of 
the purity decay is calculated in the same way. The only 
difference is that we have to consider, in addition, cross 
terms of the type ai^z ® <^i.z, i-e. terms that contain 
operators of different qubits. The contribution of these 
terms, however, vanishes when performing the trace over 
the bath variables in Eq. (|34|l . Thus, we can conclude 
that within this analytical estimate, the purity decay rate 
is identical for both the individual bath model and the 
common bath model. 

The so-called quantum degree 

Q= max (*mc|Pout |*mo) (40) 

Pout:|*mo> 

is the overlap of the state obtained after the simulated 
gate operation and the maximally entangled Bell states. 
Finally the entanglement capability C is defined as the 
smallest eigenvalue of the density matrix resulting from 
transposing the partial density matrix of one qubit. As 



8 




FIG. 6: (Colour online) Temperature dependence of the deviation of the four gate quality factors from their ideal values for the 
CNOT gate. The decoherence due to phonons is taken into account. The black line shows the results for a single phonon bath 
and the red line is for two phononic baths. The characteristic qubit energies are Ea = 1/4 GHz and the tunnel amplitudes are 
Ai = Es {i = 1, 2) due to the spacing of the double dots. In the curves for the deviation of the purity, we included lines for the 
analytical expressions 1 from Eq. I34II and 2 from Eq. 1391 . 



shown in Ref. the non-negativity of this smallest 

eigenvalue is a necessary condition for the separability of 
the density matrix into two unentangled systems. The 
entanglement capability approaches —0.5 for the ideal 
CNOT gate. 

It has been shown that the controUed-NOT (cnot) 
gate together with single-qubit operations is sufficient 
for universal quantum computation. Here, we investi- 
gate the decoherence during a CNOT gate which generates 
maximally entangled Bell states from unentangled input 
states. In FiguresElandElthe simulated gate evolution in 
the presence of phonon baths is shown. Using the system 
Hamiltonian, the CNOT gate can be implemented through 
the following sequence of elementary quantum gates^SiSS 

UcNOT = U^^ cxp (^-i^&z.ij exp (^-i^az,2j x 

X exp (^~i^azs^z,2^ exp (^-i^az,ij u'^\ 

(41) 

(2) 

where U^i denotes the Hadamard gate operation per- 
formed on the second qubit. This gate sequence just 
involves one two-qubit operation at step three. The pa- 
rameters for the numerical calculations are given below 



Figs. El and □ 

In Fig. El the gate quality factors for the case of a single 
or two distinct phononic baths are shown. It is observed 
that for the case of a single phonon bath they achieve 
better values. This offset is due to the larger number 
of non-vanishing matrix elements in the coupling of the 
noise to the spin components for the two bath case. Here, 
due to several non-commuting terms in the coupling to 
the bath and the different Hamiltonians needed to per- 
form the individual steps of the quantum gate, the gate 
quality factors saturate when the temperature T is de- 
creased. This happens at around T = = 12 mK corre- 
sponding to Eg = 1/4 GHz as the characteristic energy 
scale. 

Figure [3 depicts the same behaviour of the gate qual- 
ity factors as in FigureElwith the only difference that the 
tunnel coupling A2 is smaller by a factor of 4 during the 
Hadamard operation. The qualitative behavior is very 
similar to that in Figure O but the deviation from the 
ideal values for the gate quality factors is much smaller 
and already fulfills the criterion of an allowed deviation of 
10^^. The reduction of the tunnel amplitudes by a factor 
4 corresponds to a very small change of the distance d in 
the two qubits (namely, from 100.0 nm to 100.3 nm) ow- 
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FIG. 7: (Colour online) Temperature dependence of the deviation of the four gate quality factors from their ideal values for 
the CNOT gate. The decoherence due to phonons is taken into account. The black line shows the results for a single phonon 
bath and the red line is for two phonon baths. The characteristic qubit energies are Es — l/i GHz and the tunnel amplitude 
during the Hadamard operation on the second qubit is A2 = l/iEs, i.e., a factor 4 smaller than in Figure |S| In the curves for 
the deviation of the purity, we included lines for the analytical expressions 1 from Eq. 13411 and 2 from Eq. 1391 . 



ing to the Gaussian shape of the electron wavefunctions, 
provided their distance is sufficiently large^^. 

We have already mentioned that the phonon contri- 
bution to decoherence still allows for the fidelity values 
below the threshold 1 - F < 10"'' from Ref. 39]. For a 
reliable quantum computer, however, such intrinsic de- 
coherence mechanisms should beat the threshold at least 
by an order of magnitude. This can be achieved as fol- 
lows: As we have seen, the Hadamard gate is the step 
limiting the performance as during the Hadamard the 
system is vulnerable against spontaneous emission at a 
rate 7 oc -E^ , where E is the typical energy splitting of the 
single qubit. The duration of the Hadamard, on the other 
hand, scales as r cx 1/E. Thus, the error probability and 
the purity decay reduces to 1 — e~'>"^ ~ 7r cx i?^. Thus, 
by making the Hadamard slower, i.e., by working with 
small tunnel couplings between the dots, the gate perfor- 
mance can be increased. This works until Ohmic noise 
sources, electromagnetic noise on the gates and controls, 
takes over. This is demonstrated nicely in Fig. [7| where 
the CNOT gate for a modified Hadamard opertation (on 
the second qubit) with A2 = £2 = (l/4)i?s is depicted. It 
is clearly observed that by decreasing the tunnel matrix 
element and by increasing the evolution time the deco- 



herence is reduced and the threshold for the gate quality 
factors to allow universal quantum computation^^ can be 
achieved. 

The gate quality of a CNOT under decoherence has 
been studied in Refs. for standard collective 

and/or single-qubit noise in Ohmic environments. The 
single-qubit case for charge qubits in GaAs has been stud- 
ied in Ref. |3 with emphasis on non-Markovian effects. 
Even in view of this, and in view of the emphasis of the 
strong tunneling regime, that work arrives at the related 
conclusion that intrinsic phonon decoherence in this sys- 
tem can be limited. Please note, that the approximations 
in the microscopic model give an upper bound of valid- 
ity for the validity of effective H amilt onians as studied 
in Ref. [l^ as descibed in Refs. I7lll2lll2^ . The work 
presented here is not affected by this restriction due to 
the emphasis of the case of small tunnel coupling. 

V. CONCLUSIONS 

We have analyzed the influence of a phononic environ- 
ment on four coupled quantum dots which represent two 
charge qubits. The effective error model resulting from 
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the microscopic Hamiltonian does not belong to the fa- 
mihar classes of local or collective decoherence. It con- 
tains a dipolar and quadrupolar contribution with super- 
ohmic spectra at low frequencies, uj^ and uj^ respectively. 
The resulting decoherence is an intrinsic limitation of any 
gate performance. In particular, we have investigated 
within a Bloch-Redfield theory the relevant rates and 
the quality of a CNOT gate operation. The two employed 
models of coupling the qubits to individual heat baths 
versus a common heat bath, respectively, yield quantita- 
tive differences for the gate qualifiers. Still the qualitative 
behavior is the same for both cases. 

Within an analytical estimate for the purity loss, we 
have found that the decoherence plays its role mainly 
during the stage of the Hadamard operation. The physics 
behind this is that during all the other stages, the bath 
couples to the qubits via a good quantum number. Con- 
sequently, during these stages, the decoherence rates are 
dominated by the spectral density of the bath in the 
limit of zero frequency which for the present case of a 
super-ohmic bath vanishes. The results of our analyti- 
cal estimate compare favorably with the results from a 
numerical propagation. 

The fact that on the one hand, the bath spectrum is 
super-ohmic, while on the other hand, the Hadamard op- 
eration is the part that is most sensitive to decoherence, 
suggests to slow down the Hadamard operation by us- 
ing a rather small tunnel coupling. Then, decoherence is 
reduced by a factor that is larger than the extension of 
the operation time. This finally results for the complete 
gate operation in a reduced coherence loss. Thus, the 
gate quality is significantly improved for dots with weak 
tunnel coupling and can intrinsically meet the threshold 
for quantum error correction. 
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APPENDIX A: AVERAGE OVER ALL PURE 
STATES 

In this appendix, we derive formulas for the evalu- 
ation of expressions of the type trip A) and tr (pApB) 
in an ensemble average over all pure states p = \'ip){'4>\. 
The state is an element of an A^-dimensional Hilbert 
space. Decomposed into an arbitrary orthonormal basis 
set {\n)}n=i...N , it reads 



bution 



c„ n 



(Al) 



where the only restriction imposed on the coefficients c„ 
is the normalization = J^n 1*^"!^ ^ ^- Hence the 

ensemble of pure states is fully described by the distri- 



P(ci,...,c^) = 7jv5(l-5]|c„n. 



(A2) 



We emphasize that P(ci, . . . , cat) is invariant under uni- 
tary transformations of the state lip). The prefactor 77V 
is determined by the normalization 



y d^ci . . .(fcN P{ci, . . . ,cjv) = 1 



(A3) 



of the distribution, where j <Pc denotes integration over 
the real and the imaginary part of c. 

The computation of the ensemble averages of the coef- 
ficients with the distribution ljA2|l is straightforward and 
yields 



N 



N{N + 1) 



(A4) 
)• (A5) 



Using these expressions, we consequently find for the en- 
semble averages of the expressions tT{pA) and tv^pApB) 
the results 



tripA) = {^j\A\i;) 



tiA 
1^' 



(A6) 



tiipApB) ^ {i^\A\^){ij\B\ib) ^ 



tr(A)tr(B) + tr{AB) 



N{N + 1) 



(A7) 



which have been used for deriving the purity decay H33|) 
from Eq. 

While this averaging procedure is very convenient for 
analytical calculations, the numerical propagation can be 
performed with only a finite set of initial states. In the 
present case, the averages are computed with the set of 
36 states given after Eq. In the present case, we 

have justified numerically that both averaging procedures 
yield the same results. Thus, it is interesting whether this 
correspondence is exact. 

For the case of one qubit, N — 2, the discrete set of 
states is given by the states — ci|l) -I- C2|2) where 
(ci, C2) is chosen from the set of 6 vectors 



1 

'71 



J4> 



(A8) 



where (j) = 0, 7r/2, tt, 37r/2. Computing the averages for 
the states (IA8|I is now staightforward and shows that 
this discrete sample also fulfills the relations ljA4|l and 
(|A5|) . Thus, we can conclude that for the computation of 
averages, both the discrete and the continuous sample. 

For more than one qubit, however, arises a difference: 
While the sample of all pure states also contains entan- 
gled states, these are by construction excluded from set 
of direct products of the 6 one-qubit states ljA8|l . Still 
our numerical results indicate that the different samples 
practically result in the same averages. 
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